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ABSTRACT 

Qh' Measuring solar-like oscillations in an ensemble of stars in a cluster, holds promise for testing 

Q , stellar structure and evolution more stringently than just fitting pa rameters to sing le field stars. 

^ ■ The most ambitious attempt to pursue these prospects was bv iGillilandetall (Il993h who 

I targeted 1 1 turn-off stars in the open cluster M67 (NGC 2682), but the oscillation amplitudes 

. were too small (< 20 /imag) to obtain unambiguous detections. Like Gilliland et al.. C1993) we 

also aim at detecting solar-like oscillations in M67, but we target red giant stars with expected 
amplitudes in the range 50-500 /imag and periods of 1 to 8 hours. We analyse our recently 
. published photometry measurements, obtained during a six-week multisite campaign using 

' nine telescopes around the world. The observations are compared with simulations and with 

5^ \ estimated properties of the stellar oscillations. Noise levels in the Fourier spectra as low as 

27/xmag are obtained for single sites, while the combined data reach 19 /imag, making this 
the best photometric time series of an ensemble of red giant stars. These data enable us to 
make the first test of the scaling relations (used to estimate frequency and amplitude) with an 
homogeneous ensemble of stars. The detected excess power is consistent with the expected 
signal from stellar oscillations, both in terms of its frequency range and amplitude. However, 
our results are limited by apparent high levels of non-white noise, which cannot be clearly 
separated from the stellar signal. 

Key words: stars: red giants - stars: oscillations - stars variables: other - open clusters and 
associations: individual: M67 (NGC 2682) - Techniques: photometric. 



1 INTRODUCTION 

The firs t clear detection of so l ar-like oscill ations in a red giant star 
(^ Hya; iFrandsen et"al]|2002l : ISteIldl20oA opened up a new part 
of the Hertzsprung-Russell diagram to be explored with astero- 
seismic techniques. Follow ing that discovery, detailed analysis of 
fHya has been performed jTeixeira et al]|2003l : ISteIlo et alj|2004 
l200d) and new discoveries of oscillations in similar stars have 
emerg ed (e Oph and rj Ser; iBarban et al.l |2004| : iDe Ridder et al.l 
l2006h . These results are all based on radial velocity measurements 



of high precision (cr ~ 2 m/s) but from non-continuo us observa- 
tions, which imposes large a mbiguities on the results l lStello et"ail 
2006; iDe Ridder etal] l2006h . With oscillation periods of a few 
hours, these stars require a time base of roughly one month, which 
can only be obtained on small telescopes. But the current lack of 
high-precision spectrographs on small telescopes makes a multisite 
campaign impossible. 

However, using photometry makes it feasible to incorporate 
many l-2m class telescopes in a multisite campaign, and it fur- 
thermore provides the possibility to observe many stars, like in 
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(B-V) 



Figu re 1. Colour-magnitude d iagram of the open cluster M67 (photometry 
from lMontgomerv et alj 19931) . The target stars are indicated with identifier 
numbers con'esponding to those given in Tables [T] and |2] The solid line is 
an isochrone with (m - M) = 9.7 mag. Age = 4.0 Gyr Z = 0.0198 
and Y = 0.2734 from the BaSTI database jPietrinfemi et alj2004 . 



a cluster, simultaneously. Detecting oscillations in a set of clus- 
ter stars potentially increases the power of the asteroseismic mea- 
surements due to the additional constraints provided by the com- 
mon parameters of the cluster members (age and composition). 
Until recently, there has been no such ground-based photomet- 
ric campaign aimed at detecting solar-like oscillations in red gi- 
ant stars. A number of attempts have been made to detect oscil- 
lations in more Sun-like stars (hotter and less lu minous) of the 
open cluster M67 (Gilliland & Bro wn 1988; Gillil and etalj|l99ll : 
iGilliland & BrowrJll992h, with the most ambitious multisite ef- 
fort made bv lOilliland et all ( Il993h . Despite noise levels as low 
as 0.29 mmag per minute integr ation no una mbiguous detections 
were claimed. In a recent paper dStello et alj|20 06 hereafter Paper 
I), we reported observations from a large six-week multisite cam- 
paign also aimed at M67. However, unlike the previous studies our 
campaign was optimized for the slightly brighter and longer period 
red giant stars (see Fig. and also covered a much longer time 
span. The long oscillation periods mean that non-white noise such 
as drift is more crucial for this project than in the previous studies. 
The data set reported in Paper I was based on very different sites, 
many with unknown long-term stability performances. A realistic 
estimate of the final non-white noise in the data could therefore not 
be obtained prior to observations. 

In this paper our main emphasis is on the time series analysis 
of the red giant stars (Sect.|4ll, based on the data described in Paper 
I. We report in Sect. |2] on an additional independent data reduc- 
tion method to further obtain lower noise in the Fourier spectra. In 
Sect. [3] we estimate the oscillation characteristics and simulate in 
Sect. |5] the expected outcome for each target without the presence 
of non-white noise to facilitate the analysis of the observations. We 
give our conclusions in Sect.|6] 



2 OBSERVATIONS AND DATA REDUCTION 

The data are from a global multisite observing campaign with nine 
0.6-m to 2.1-m class telescopes from 6 January to 17 February 
2004 (Paper I). The photometric time series of those stars within 
the field-of-view of all telescopes comprises roughly 18000 data 
points. 

After calibrating the CC D images we used the momf package 
jKieldsen & Frandsedl 19921) . which calculated differential time se- 
ries of 20 red giants relative to a large ensemble of stars (from 116 
to 358 stars, depending on telescope field-of-view). We performed 
the following three initial steps to improve the signal-to-noise in 
the final Fourier spectra of the time series: (I) sigma clipping, (2) 
correction for colour extinction and (3) calculation of weights for 
each data point. For further details see Paper I. 

In addition to the time series produced by D.S., as described 
above, H. B. construc t ed dif ferential time series following the ap- 
proach of lHonevcutj h992l) . This was based on the raw time se- 
ries (also calculated by momf), but using only the red giant stars 
as the ensemble. For each target star, the differential time series 
was calculated by subtracting a reference time series that did not 
include the star itself. The reference time series comprised offsets 
for each data point (CCD image), which were a weighted average 
of the ensemble. For each data point the weight w given to each 
star was calculated as w = l/(crptp + cTmin), where dptp is the 
local point-to-point scatter and (Tmin = 1 mmag is a fixed mini- 
mum noise value to prevent a single star with very low noise from 
dominating the ensemble. Using the relatively small homogeneous 
ensemble comprising only red giants has the advantage of better re- 
moving colour-dependent extinction, and hence providing a lower 
noise level in the Fourier spectrum at low frequencies. However, 
simulations showed that the reference time series will include stel- 
lar oscillations with amplitudes up to 30% of those we want to de- 
tect. The stars with the longest periods (hence largest amplitudes) 
will be the most affected because the sample is dominated by these. 

In most stars, the two methods produced very similar noise 
levels in the time series, within 10%. However, on a few stars dif- 
ferences of up to 50% were seen. In the following, we chose the 
time series for each star and each site with the lowest noise in the 
Fourier spectrum in the interval 300-900 /iHz, which is outside the 
frequency range where the stars are expected to oscillate. 



3 EXPECTED OSCILLATION SIGNAL 

To better interpret our results, we have estimated the characteris- 
tics of solar-like oscillations expected in the red giant stars. We 
scaled the oscillation parameters of the Sun to predict amplitude, 
central frequency of excess power and the large frequency separa- 
tion. These predictions are used in Sects. |4]and|5]to compare with 
the observations. 

The predicted amplitude in the Johnson V filter 
(Acen = 544 nm) was derived using the scaling relation by 
lKieldsen~& Bedding ( 1995J (using 1 ppm = 1.086 /imag): 

L/Lq(5.1 ±0.3) /xmag 



{5L/L)^ = 



(A/550 nm) (T^ff /5777 K)2 (M/Mq ) 



(1) 



These amplitudes were used as a guide while planing the observa- 
tions. However, recent theoretical studies indicate that the L/M- 
scaling may over-estimate the amplitude for main-sequence stars 
and that (L/M)°'^-scali ng might provide a more realistic predic- 
tion jSamadi et aLll2005h . which we will take into account in eval- 
uating of our results. We note that extrapolating amplitudes for red 
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Table 1. Properties of red giant target stars. 



No 




B~V 


L/Lq 




SL/L /imag^ 


^max 


Al/Q 


Cross-ref^ 


pd 




mag 


mag 




K 


(L/M) 




^Hz 


/jHz 




% 


3 


9.72 


1.37 


250.6 


3920 


2080 


434 


3.7 


0.7 


S978/G8 


53/95 


11 


9.69 


1.36 


243.5 


3960 


1980 


417 


3.9 


0.8 


S1250/G4 


58/95 


4 


10.30 


1.26 


87.0 


4330 


592 


170 


15.1 


2.2 


S1016/G5 


78/93 


21 


10.47 


1.12 


51.9 


4727 


296 


99 


34.4 


4.2 


SI 592/ 


48/92 


8 


10.48 


1.11 


50.8 


4750 


287 


97 


35.8 


4.3 


S1010/G2 


82/96 


9 


10.48 


1.10 


50.2 


4772 


281 


95 


36.8 


4.4 


S1084/ 


80/92 


10 


10.55 


1.12 


48.2 


4727 


275 


94 


37.0 


4.4 


S1279/G7 


79/92 


20 


10.55 


1.10 


47.1 


4772 


264 


91 


39.2 


4.6 


S1479/ 


76/95 


2 


10.59 


1.12 


46.4 


4727 


265 


92 


38.4 


4.6 


S1074/ 


74/91 


18 


10.58 


1.10 


45.8 


4772 


256 


89 


40.3 


4.7 


S1316/ 


73/95 


16 


10.76 


1.13 


40.3 


4703 


232 


84 


43.5 


5.0 


S1221/ 


92/90 


5 


11.20 


1.08 


25.4 


4815 


140 


58 


74.8 


7.6 


S1054/G9 


93/64 


17 


11.33 


1.07 


22.4 


4835 


122 


53 


86.0 


8.4 


SI 288/ 


94/96 


7 


11.44 


1.06 


20.2 


4854 


109 


49 


96.9 


9.2 


S989/G12 


95/95 


19 


11.52 


1.05 


18.7 


4873 


101 


46 


106.0 


9.9 


S1254/ 


94/95 


15 


11.63 


1.05 


16.9 


4873 


91 


43 


117.3 


10.6 


S1277/ 


95/95 


14 


12.09 


1.01 


11.2 


4945 


58 


31 


187.3 


15.2 


S1293/ 


96/93 


12 


12.11 


1.01 


11.0 


4947 


57 


31 


190.2 


15.4 


S1264a/G15 


0/75 


13 


12.23 


0.99 


9.9 


4966 


51 


28 


213.2 


16.8 


SI 305/ 


96/95 


6 


12.31 


0.99 


9.2 


4971 


48 


27 


230.2 


17.8 


SI 103/ 


95/— 



Photometry from lMontgomerv at aljfl993l) . 
^ Estimated amplitudes based on both [LIM") a nd {L/My^'"^ scaling. Note t hat 1 ppm = 1.086 /imag. 
IDs starting with an S are from lSandersI i\91% . while G are from Gilliland etalJ 1199 ih . 

Mem bership probabilities are from lZhao et"ai] h993l) and lSandersi |1977|) respectively. All targets have high probabilities 
inloirard et al. 1 1989) (P > 95%). 



giant stars based on these scaling relations is very uncertain, and 
has so far not been thoroughly tested by calculations of theoretical 
pulsation models or observations of these stars. 

The characteristic frequency domain within which a star is os- 
cillating was estimated as the central frequency of the excess power, 
which was obtained by scaling the acoustic cut-off frequency of the 
Sun cBrown etal.199 1) 



M/Mq 



{R/RqY Vr,ff/5777K 
M/Mq (reff/5777K)^-5 



X 3050 AiHz 



X 3050 ^Hz , 



(2) 



using L/Lq = (7?/R0)2(Tcft/5777K)''. This scaling relation 
gives very good agreement with the frequency range of the solar- 
like oscillations observed in main sequence stars and also in red 
giants Igedding & Kjeldsen 2Q^. 

Finally, we predict the expected frequency spacing, Afo, 
between modes of the s ame degree in the power spectrum 
jKieldsen & Beddindll995h 

^''o = \ J,^ ,\ ^ X (134.92 ± 0.02) /iHz (3) 



(i?/R0)i-5 

(M/M0)°-^(reff/5777K)^ 
(L/LQ)o.r5 



X (134.92 ± 0.02) /iHz , 



To calculate these parameters we made rough estimates of 
the stellar mass, luminosity, and effective temperature, which 
are summarized in Table [T] For all stars we adopted a mass 
of Af — 1.35M0 corresponding to the mass at the base of 
the red giant branch for an isochrone with Age = 4.0 Gyr, 
Z = 0.0198 and Y = 0.2734 (BaSTI database: Ipietrinfemi et al] 
which matches the cluster colour-magnitude diagram (see 



Fig.[T). This is also in good agreement with the turn-off mass by 
VandenBerg & Stetson (2004). We derived L and T^s using {V, 
B — l^)-photometry of the cluster (Montgome ry et alll 19931) and 
interpolation of the BaSeL grid dLeieune et aL, 19981) . For our pur- 
pose it was sufficien t to adopt a typical surface gravity for all red 
giant s of log g=2.5 ( lAllenlll973h , [Fe/H]=0.0 in agreement with 



,Nissen et al 



and (m — M) = 9.7 mag corresponding 
to d = 870 pc, which is within 5% of previous investig ations 
(e.g. iMontgomerv et alJ[T993l : IVandenBerg & Stetsorj|2004l) . The 
assumption of a common log g affected our tempareture estimates 
by less than 100 K and our luminosity estimates by up to 10% for 
the stars investigated. We note that, because the temperatures of 
our target stars are roughly the same, the expected amplitudes are 
proportional to 1 / Vmax 



V5777K 



550 nm 
A 



(10.1 ±0.8) mmag.(4) 



4 TIME SERIES ANALYSIS 

In Table 12 we give for each star the noise levels (in amplitude) in 
the Fourier spectra measured in two frequency intervals based on 
single-site data. The noise denoted (tiooo-3000 hHz represent the 
lowest noise level (white noise), while asoo-goOfiHz is the noise 
level closer to, but still outside the expected frequency range of the 
oscillations. The best data were from SAAO, La Silla and Kitt Peak, 
with mean noise levels of roughly 40 /imag for the best stars. How- 
ever, at SAAO only a few stars were observed due to the small field- 
of-view. The Kitt Peak data generally had a slightly lower point- 
to-point scatter but were more affected by extinction and hence 
showed more drift noise. As we seek to optimize the signal-to-noise 
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Table 2. Mean noise level (in /jmag) in the Fourier spectrum for two frequency intervals (stars plotted in Fig.|2]are in boldface). Site abbre- 
viations are: SSOi (Wide Field Imager at Siding Spring Observatory, Australia); SSO2 (Imager at Siding Spring); SOAO (Sobaeksan Optical 
Astronomy Observatory, Korea); SAAO (South Africa Astronomical Observatory); RCC (Ritchey-Chretien-Coude at Piszkesteto, Konkoly Ob- 
servatory, Hungary); Sch (Schmidt at Piszkesteto); LaS (La Silla Observatory, Chile); LOAO (Mt. Lemmon Optical Astronomy Observatory, 
Arizona); Kitt (Kitt Peak National Observatory, Arizona); Lag (Mt. Laguna Observatory, CaHfornia). 



No SSOi SSO2 SOAO SAAO RCC Sch LaS LOAO Kitt Lag 

O" 1000-3000 mHz 1 O" 300-900 mHz 



3 


291, 


372 


298. 311 


400, 535 


-, 


- 


-, - 


809, 1090 


481, 585 


349, 545 


-, 


- 


604, 746 


11 


2930, 3325 


231, 235 


334, 366 


294, 


417 




1096, 1737 


102. 


120 


476, 688 






152, 196 


4 


99, 


110 


138, 134 


147, 137 


43, 


48 


248, 386 


139, 


144 


39, 


41 


164, 268 


62, 


91 


123, 173 


21 






229, 260 


359, 394 








227, 


302 






149, 153 








8 


101, 


87 


155, 155 


117, 134 


48, 


52 


171, 295 


117, 


111 


37, 


38 


146, 229 


47, 


59 


118, 152 


9 


94, 


95 


154, 174 


129, 141 








120, 


130 


40, 


42 


194,315 


27, 


31 


118, 157 


10 


82, 


69 


97, 122 


119, 115 


68, 


78 


164, 284 


113, 


124 


37, 


40 


147,217 


34, 


42 


93, 124 


20 


86, 


83 


115, 119 


152, 157 








150, 


196 


39, 


39 


258, 465 






191,336 


2 


96, 


101 


135, 157 


148, 151 








138, 


138 


53, 


52 


150, 202 


32, 


38 


125, 166 


18 


77, 


78 


122, 134 


187, 176 








120, 


136 


33, 


36 


247,451 






156, 280 


16 


92, 


106 


167, 164 


259, 239 








150, 


150 


46, 


48 


289,312 






184, 295 


5 


96, 


87 


160, 171 


213, 222 


45, 


47 


219, 324 


165, 


186 


57, 


59 


108, 128 


43, 


51 


114, 175 


17 


100, 


103 


126, 148 


199, 194 






255, 341 


162, 


168 


101, 


122 


154, 303 


36, 


43 


163, 225 


7 


108. 


95 


159, 178 


250, 235 


66, 


81 


319, 518 


171, 


170 


46, 


50 


122, 160 


39, 


47 


511, 810 


19 


113. 


103 


160, 180 


219, 202 








182, 


180 


55, 


64 


142, 235 


138, 


230 


153, 191 


15 


109. 


91 


153, 145 


237, 273 






331,453 


175, 


225 


163, 198 


101, 116 


45, 


53 


131, 184 


14 


109. 


94 


178, 192 


330, 340 






414, 563 


205, 


237 


50, 


51 


117, 129 


63, 


67 


131, 185 


12 


231. 


240 


486, 567 


719, 692 


140, 


156 


471,470 


198, 


227 


276, 317 


371,587 






224, 299 


13 


111. 


104 


187, 192 


349, 336 








224, 


227 


111, 128 


96, 119 


74, 


85 


148, 196 


6 






231,260 


351, 378 








232, 


289 






161, 186 









in the final Fourier spectra (in amplitude) based on the combined 
data from all sites, the data from La Silla and Kitt Peak will domi- 
nate for the majority of the stars. For further details about the noise 
properties see Paper I. 

Due to nightly drifts in the data we saw very strong peaks 
at l-^ cycles per day (corresponding to 11.57, 23.15, 34.72 and 
46.30 /iHz) in the Fourier spectra based on individual sites. Even 
when combined, the data still showed significant excess power due 
to these drifts. This was a serious problem because the most promis- 
ing stars in the ensemble are expected to oscillate in the affected 
frequency range. We decided to remove ("clean") these specific 
frequencies using standard iterative sine-wave fitting, on a site-by- 
site basis. Compared to a classic high-pass filter with a smoothly 
varying response function, this method has the advantage that it re- 
moves only a small amount of power and only in a very limited 
and well-defined frequency range, which was important because of 
the expected low frequencies of the oscillations in the stars. The 
mean noise levels in the Fourier spectra at 300-900 /iHz were re- 
duced by 2-13% (in amplitude) as a result of this cleaning process. 
We did not decorrelate the time series against external parameters 
(e.g. airmass, sky background, position on CCD) because we found 
that to have too dramatic and uncontrolled effects on the time series 
on time scales similar to the expected stellar oscillations. This was 
based on decorrelation of simulated time series. 

After the initial cleaning of the dominant low frequency noise 
peaks, we calculated the Fourier spectra of the combined data to 
search for excess power. The Fourier spectra were divided by the re- 
sponse function of the cleaning process to restore the overall power 
distribution. The response function was obtained by performing the 
same cleaning process on simulated time series of white noise with 
the same window function as the observations. The final response 
function for each star and each site was the average Fourier spec- 



trum based on 100 simulations. In Fig. [2] we show the spectra of 
the best twelve stars (noise below 50 /imag in the frequency range 
300-900 ^Mz). Three other stars (Nos. 16, 19 and 20) fulfill this cri- 
terion but have been omitted due to significantly higher noise levels 
compared to stars of similar brightness. The solid white line in each 
panel of Fig.|2]is the smoothed spectrum, which was obtained by 
smoothing twice, with a boxcar width of 30 /^Hz followed by one 
of 10 /iHz. The frequencies where we expect the stellar oscillations 
(i^max, Eq. I2J are indicated with the white downward-pointing ar- 
row head, and the dotted horizontal line shows (8(7300-900 mHz)^. 
The increase in noise towards low frequencies, are due to drifts 
in the data from instrumental and atmospheric instabilities as dis- 
cussed in Paper 1. 



4.1 Location of excess power 

We expect the oscillation modes to be located at lower frequen- 
cies for the more luminous stars (see Eq.[2]and Table[Tll. To search 
for general trends in the Fourier spectra, we grouped the stars ac- 
cording to luminosity. We formed three groups: the clump stars (5 
stars), the stars between the clump and the lower RGB (4 stars), 
and the lower RGB stars (2 stars), but we excluded star No. 4 due 
to its sole location in the colour-magnitude diagram (see Fig. [T]!. 
For each group we averaged their Fourier spectra and smoothed 
the final spectrum. Smoothing was done twice. First, with a wide 
boxcar (width=100 /iHz) to smear out humps of power originating 
from the individual spectra to better illustrate the overall distribu- 
tion of power within each group. We then used a second boxcar 
(width=10 ^Hz) to smooth small point-to-point variations in the fi- 
nal plot. The result, shown in Fig. [3] confirms the general trend 
of shifting power excess as a function of luminosity, in agreement 
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Frequency/^Hz 
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Figure 2. Fourier spectra of red giant stars (identifier shown in each panel). The solid white line is the smoothed spectrum, and the white arrow head indicates 
where we expect the stellar excess power (Table[TJ. The dotted horizontal line shows (3(T3oo-900 ;tHz)^- The inset shows the spectral window on the same 
frequency scale as the main panel. 



with expectations (see arrows in Fig. [5] showing expected location 
of excess power). 

We investigated whether the observed power shift could have 
been produced by the reduction methods, (i) We inspected Fourier 
spectra and weights from each site and concluded that the power 
shift was not an artifact originating from high weights given to a 
few sites, (ii) We reanalysed the data without removing the peaks 
at 11.57, 23.15, 34.72 and 46.30 /iHz, but saw the same trend, al- 
though with more power at low frequencies, (iii) Finally, we in- 
vestigated whether our data reduction could introduce such a trend 
due to more pronounced drift noise for brighter and cooler stars. 
The best comparison stars for this purpose were the blue stragglers, 
which had similar brightness to the red giants and hence similar 
inp ut parameters for the extrac tion of the photometry (for details 
see iKield sen & Frandsenll993 and Paper I). However, their colours 
were very different from the red giants and from star to star, and 
only a few blue stragglers were observed at most sites. The blue 



stragglers showed different levels of excess power at low frequen- 
cies from star to star, but no causal relation could be found. 

In summary, we found no evidence for a non-stellar source 
that could explain the observed pattern of excess power. The ob- 
served behaviour is exactly as expected for oscillations. However, 
due to the lack of a good ensemble of reference stars, similar to 
the red giant sample, we could not exclude with certainty that the 
photometric data reduction might have caused the shift as a result 
of more drift noise for more luminous and cooler stars. 



4.2 Amplitude of excess power 

To search for further evidence of solar-like oscillations, we mea- 
sured the amount of excess power above the background noise, and 
converted that into amplitude per mode in order to compare it with 
expected amplitudes. 
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Figure 3. Average power distribution as a function of frequency for three 
groups of stars. The arrows show the mean position of the expected excess 
power, (fc'max), for the stars (stellar IDs in brackets) that have been used to 
calculate the corresponding power distribution (see text for details). 



The noise level in the power spectrum, as a function of 
frequency, piy), was estimated using the following expression: 
p{v) — ajv^'"' +Pwn- We tried v to the power of —1, —1.5 and 
—2, and found —1.5 provided the best fit in the region around the 
excess power We measured the white noise, pwn, as the mean level 
in the frequency range 300-900 pHz. Then, having fixed Pwn, we 
calculated a requiring that the integral of piy) in the range 50- 
900 /iHz should equal that of the power spectrum. Assuming the 
excess power above this noise fit is real, we estimated th e ampli - 
tude per mode based on the approach of iKieldsen et al] dlOOSh . 
These amplitude estimates are independent of the mode lifetime, 
but require an assumption of the number of modes that are excited. 
We have estimated the amplitudes for two extreme scenarios: (1) 
only radial modes are excited and, (2) three additional non-radial 
modes per radial mode, neglecting any difference in mode visibili- 
ties due to varying cancellations in full-disk observations. The ap- 
proach is illustrated in Fig. |4] for three stars that all show a hump 
of power close to the expected frequency for solar-like oscillations 
(see arrow in left panels). The power spectra (Fig.[2j were divided 
by the integrated power of the spectral window to obtain power 
density spectra (Fig.|4]left panels). The background noise was sub- 
tracted from the smoothed spectra and the resulting power density 
was multiplied by the mean mode spacing (Aj/q for only radial 
modes, and Aj/q /4 for both radial and non-radial modes). We then 
took the square root to convert to amplitude (Fig.|4]right panels). 

The amplitude at i/max for each star is plotted versus 
L/{MT'^f^) in Fig. [5] and compared with results from simulations 
which used input amplitudes derived from Eq. [T] The filled dots 
corresponds to the values read from the solid black curves in Fig.|4] 
(right panels), while the empty circles in Fig.[5]corresponds to the 
solid grey curves (Fig. |4). However, for stars expected to oscillate 
at frequencies lower than ~ 50 /iHz, we could not make realistic 
estimates of the noise contribution due to the apparent rapid rise in 
noise at low frequency, making it impossible to disentangle noise 
and stellar signal. We could only obtain upper limits, which are in- 
dicated with arrows in Fig. [5] Since these are upper limits, we only 
show the value derived assuming only radial modes are excited in 
these stars. 

Calculations by IChristensen- Dalsgaard | |2004|) and 
iDziembowski et al] ( 12001 ) on more massive and luminous 
stars than our targets suggest that only radial modes are excited 



Figure 5. Measured amplitude versus L/(MT2jj) for the best target stars 
(stellar IDs are indicated). Arrows show observed upper limits for stars 
where we could not estimate the noise contribution. For the other stars the 
filled dots (assuming only radial modes) and empty circles (assuming ad- 
ditional three non-radial modes) bracket the inferred amplitudes. The solid 
black line shows the mean amplitude from 100 simulations using Eq. [T] 
as input amplitude. The dotted lines are ±2(7, where cr is the rms scatter 
in the simulations. The grey line indicates the mean amplitude based on 
(L/A/)f'7-scaling 



to observable amplitudes in red giant stars. No similar theoretical 
pulsation analysis have been published for the red giants in M67 in 
order to determine whether non-radial modes are expected. How- 
ever, stars Nos. 13 and 14 are expected to oscillate at frequencies 
that are significantly higher than for typical red giant stars, and 
are sim ilar to those of sub giants. Obs ervations of the sub giants 
7? Boo IKieldsen et al. 2003), /3Hyi ( iBedding et d] |2007|) and 
1/ Ind ( Carrier et al.„2007„) clearly show evidence of both radial and 
non-radial modes. 

The results in Fig. |5] show that most stars are in agreement 
with the scaling relation in Eq. [T] (illustrated with the black line). 
There is a tendency that the scaling based on (L/Af (grey line) 
is in less good agreement with the data. However, we are not able 
to exclude either due to the uncertainty whether non-radial modes 
are present, and further because we only have upper limits for the 
most luminous stars. We note that for star No. 14 some non-radial 
modes are required to make the measured power excess in agree- 
ment (within 2a) with the scaling relation in Eq.[T] We further note 
that star No. 13 fall above the 2a region even in the extreme case 
of three additional non-radial modes per radial mode. Hence, if the 
observed hump of power for this star is due to stellar oscillations, 
it is likely that non-radial modes are excited, and that we observed 
the star in a rare state of high amplitudes. The reason for this high 
excess power could also be due to more complicated drift noise in 
the data than described by our noise fit. 

In summary, we see good evidence for oscillation power in 
three out of six stars that are expected to oscillate with frequencies 
larger than 50 /iHz. For stars expected to oscillated below 50 ^Hz, 
we were only able to provided upper limits due to difficulties in 
determining the noise levels in their power spectra. In general we 
see agreement with the predicted amplitudes from L/Af-scaling. 



4.3 Autocorrelation of excess power 

To search for a regular series of peaks, which is a typical signa- 
ture of solar-like oscillations, we calculated the autocorrelation of 
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Figure 4. Left panels: Power density spectra of three selected red giant stars. Tlie thick white solid Une is the smoothed spectrum and the thin grey line is 
a fit to the noise. The aiTow indicates the expected value of I'max- Power density was calculated from Fig.|2]by dividing with the integrated power of the 
spectral window. Right panels: Estimated amplitudes for each of the three stars, derived by subtracting the noise from the smoothed power density spectrum, 
multiplying with Auo/n (n = 1 for radial modes only, n = 4 for radial plus non-radial modes) and taking the square root. 



all the Stars in the frequency range where we expect stellar oscilla- 
tions. For most stars the autocorrelation is dominated by the aliases 
at Ic/d. In about half the stars there is a hint of a peak that coin- 
cides with the expected large frequency separation (Eq. |4j and of 
those, there are two (Nos. 7 and 13) that show intriguing results. In 
the former, the most prominent peak is located at 8.8 /xHz, which 
is very close to the expected value (Table [T). In Fig. [6] we show 
the autocorrelation of star No. 13. The autocorrelation was calcu- 
lated in the frequency range 135-285 /jHz setting all values in the 
Fourier spectrum lower than a threshold of 20-300-900 mHz (in am- 
plitude) equal to the threshold value. The most significant peak that 
is not a daily alias is located at 19.3 /xHz, which is within 15% of 
the expected value (Table [T]i. 



4.4 Stellar granulation 

In order to estimate the signal from background granulation we use 
a simple scaling of the solar granulation background observed by 
SoHO ( VIRGO green channel). The scaling uses the Harvey model 
l lHarvevHl985h to describe the size and shape of the granulation 
power spectrum. The Harvey model needs as input the characteris- 
tic time scale for granulation, as well as the granulation scatter in 
the time series. We calibrated those two parameters using the solar 
data, using the following simple assumptions: 

(1) The size of a granulation cell on the surface of a star is pro- 
portional to the scale height of the atmosphere at the stellar surface. 
The variance of the granulation in the time series is then assumed to 
be inverse proportional to the total number of cells on the surface. 

(2) The characteristic time scale of granulation is approxi- 
mated to be proportional to the ratio between the size of the individ- 
ual cells and the velocity of the cell. We estimate the granulation 
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Figure 6. Autocorrelation of star No. 13 (see text). Dashed lines indicate 
integer multiples of one cycle per day (1 1.574 ^Hz). 
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Figure 8. Autocorrelation of a simulation of star No. 13 (see text). Dashed 
lines indicate integer multiples of one cycle per day (1 1.574 /xHz). 



cell velocities to be proportional to the sound speed in the stellar 
atmosphere. 

(3) We assume that the atmosphere is isothermal (allowing a 
simple scaling of the atmospheric scale height). 

More details can be found in Kjeldsen & Bedding (in prepara- 
tion). Using these assumptions, we are able to calculate granulation 
power spec tra that agrees both in shape and level with detailed sim- 
ulations bv lSvensson & Ludwig| ( l2004l) . In Fig.|7]we plot smoothed 
power density spectra of the target stars, together with the estimated 
signal from background granulation. For the most luminous stars 
granulation is likely to contribute significantly at low frequencies, 
while for the most faint stars in our sample granulation can be ex- 
cluded as a possible cause of increasing noise. From this figure we 
can also see that there is a distinct difference in time scale between 
the granulation and the humps of power that appear to be due to 
stellar oscillations. Background granulation will therefore not pre- 
vent detection of solar-like oscillations. 



5 SIMULATIONS 

To interpret the results shown in Sect. |4] we made a series of 
simulations of ea ch target star using the method described by 
IStello et al.l ( l2004t) . We chose a regular series of input frequen- 



cies with a separation of A^o (Table [TJ. Their relative amplitudes 
were determined by a Gaussian envelope with a height correspond- 
ing to the L/A/-scaling values in Table [T] (note that (L/M)°-^- 
scaling predicts amplitudes that are roughly 0.3-0.5 times that for 
our targets). The envelope was centered at i^max (Table [T]! with a 
width equ al to 0.48 i/'max, w hich was calibrated to the observations 
of 5 Hya dStello et al.l | 2004 l). The envelope reprodu ced by this ap- 
proach is in good agreement with (3 Hyi and the Sun (Kjelds en et al.l 
I2OO5.) . The number of input frequencies was the closest integer of 
7 X 0.48 t'max/ Afo, which ensures that we have frequencies within 
the entire envelope. Because the non-white noise is difficult to esti- 
mate precisely, we chose to include only white noise in the simula- 
tions, which we set equal to the mean level in the Fourier spectra at 
300-900 ^Hz of the observed data. For each target star we ran sim- 
ulations in two parallel runs, one with a mode lifetime of 20 days, in 
agreement with the theoretical calculations on the red giant ^ Hya 
(Houdek & Gough 2002), and th e other of 2 days in agreements 
with the observations of that star jStello et al.l2006l) . 

We show here the results for the short mode lifetime (Fig.|9]l. 
The difference between the long and short mode lifetime is dis- 
cussed below. It is important to stress that the detailed characteris- 
tics of the excess power in the Fourier spectra depend very much 
on the complicated interaction between the spectral window and 
the oscillation modes (their frequencies, amplitudes and lifetime), 
as well as the random number seed. For example, using the same 
random number seed, the clump stars (Nos. 8, 9, 10, 2 and 18), 
which had different spectral windows but only slightly different 
amplitudes and input frequencies, show very different power ex- 
cesses (Fig. |9]l. For star No. 2 we show in Fig. [TO] (left panels) 
five spectra based on different random number seeds. Again, we 
see a large variation in the result, which is in agreement with Fig.|5] 
(dashed lines). This large intrinsic variation was also pointed out by 
IStello etal.l l l2004l2006h . The right panels show the corresponding 
spectra based on the long mode lifetime (=20 days). From this, we 
see that a long mode lifetime in general provides higher peaks, but 
the differences are comparable to the difference that arises from 
using a different seed. 

A comparison of observations and simulations (Figs.|2]and|9]l 
indicates that a significant part of the excess power seen in the ob- 
servations might be of stellar origin. The excess power in the obser- 
vations can be associated with either humps of power (see Nos. 13 
and 17 in Fig.[2ll or a broad "shoulder" on top of the drift noise (see 
No. 2). It is clear from the simulations that observing excess power 
in one star does not imply that we necessarily would see a similar 
hump in a similar star (or even in the same star if observed at an- 
other epoch; Fig.llOll. We note that the simulations shown in Fig.|9] 
only included radial modes. If the stars have non-radial modes these 
simulations underestimate the stellar excess power. This could be 
significant at le ast for the less evolved stars (Nos. 13 and 14). Very 
recently, iHekker et alj ( l2006h reported that non-radial modes are 
present in two red giant stars (e Oph and rj Ser) that oscillate at ap- 
proximately 60 /iHz and 130 /iHz, respectively, and hence are com- 
parable to stars in our sample. Here we also note that the scaling 
relation that was used to estimate the amplitudes is quite uncertain 
for red giant stars. 

Finally, we calculated the autocorrelation of simulated data of 
star No. 13 that reproduced the observed excess power. This was 
done by setting the input amplitude of the simulations to 200 /imag 
and including only radial modes in accordance with Fig.|4] This 
time we set the white noise equal to the value estimated at 200 /.iHz 
from our noise fit in Fig.|4] The autocorrelation was based on simu- 
lations with a short mode lifetime (= 2 days) because they produced 



Multisite campaign on M67. II. 



9 




Figure 7. Smoothed power density spectra of red giant stars (identifier shown in each panel). The dashed curve is the estimated granulation signal based on 
scaling from the Sun (see Sect l4.4t . The dotted hne indicates where stellar oscillations are expected (Eq.[2)- 



Fourier spectra more like the observed spectrum. In total nine inde- 
pendent simulations were made. In roughly half the cases we saw 
a peak in the autocorrelation of the same significance as in the ob- 
servations (Fig. |6l(. In Fig. [8] we show one of the examples which 
show a peak. However, these peaks were not at exactly the same 
frequency but appeared in an interval from approximately 16/iHz 
to 19 nHz. Some of these (as in Fig.[8]l were close to halfway be- 
tween 1 and 2 cycles per day (17.4 /xHz), and might therefore be 
caused by aliasing. In the other cases the only clear peaks were 



those clearly originating from the daily aliases at 1, 2 and 3 cycles 
per day. 

To investigate whether the peak seen at 19.3 /iHz in the auto- 
correlation of the observations was caused by interaction between 
the spectral window and random peaks we repeated the simula- 
tions using randomly spaced frequencies. For each of the three sets 
of random frequencies we made nine independent simulations, and 
inspected the autocorrelation. The differences seen between each 
of the independent simulations within each set of frequencies were 
small compared to the difference between different sets. In most 
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Figure 9. Fourier spectra of simulated red giant stars (identifier shown in each panel), with a mode lifetime of 2 days and amplitudes according to L/A/-scaling 
(Eq.[T). The solid white line is the smoothed spectrum, and the white aiTow head indicates where we expect the stellar excess power (Table[T). The inset shows 
the spectral window on the same frequency scale as the main panel. The dotted line is (30-300-900 ^iHz)^- 



cases we did not see a clear peak in the range 16-19 /iHz, but in 
a few cases a clear peak halfway between 1 and 2 cycles per day 
(17.4 ^Hz) was seen. In some cases peaks as significant as the ob- 
served 19.3 fiYlz peak were seen at other frequencies. Hence, we 
conclude that the peak at 19.3 ^Hz in the observations could be 
caused by randomly distributed peaks. 



6 CONCLUSIONS 

We have analysed high-precision photometric time series of 20 red 
giant stars in the open cluster M67. The data, recently published in 
Paper I, were from a large multisite campaign. 

We see evidence of excess power in the Fourier spectra, shift- 
ing to lower frequencies for more luminous stars, consistent with 
expectations from solar-like oscillations (see Fig. [3}. 

We further estimated the mode amplitude from the observed 



power excess assuming that it originates from solar-like oscilla- 
tions. These estimated amplitudes are in agreement with the results 
from simulations using estimated amplitudes by scaling from the 
Sun (Fig. Is}- These results further support that we have seen ev- 
idence for solar-like oscillations in some of the red giant stars in 
M67. We note that the measured amplitudes from our simulations 
show a scatter of roughly 40% from o ne epoch to the next, which 
is similar to what we see in the Sun l lKieldsen et alll200'4) . This 
intrinsic scatter should always be taken into account when evalu- 
ating observed amplitudes of solar-like oscillations. Without firm 
knowledge of whether non-radial modes are present in these stars 
and with only upper limits on the amplitudes of the most luminous 
stars, the large scatter in the measured amplitude limits our abil- 
ity to clearly pinpoint a favoured scaling relation for amplitude. 
However, we see an indication of a better match between observa- 
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tions and the L/ M-scaling iKieldsen & Bedding| | l995l) than with 
(L/Af)'''^-scahng jSamadi et alfcOOSh . 

In many stars we see apparently high levels of non-white 
noise, but its temporal variation is unknown and could not be decor- 
related. In the Fourier spectrum this noise seems to be significant at 
frequencies below 100 /iHz. This is unfortunate because the target 
stars in which we expected the most clear detections are expected 
to be oscillating in the affected frequency range. We are therefore 
not able to disentangle the noise and stellar signal in the analysis of 
these stars. Hence, the data did not allow the clear detecti ons that 
we had anticipated. Simple scaling of the solar granulation jHarvevI 
1985) shows that the power at low frequencies (lower than the ex- 
pected p-mode range) is not coming from granulation for the fainter 
stars on the low RGB, and is most likely instrumental (see Fig.|7]l. 
However, it could be granulation noise in the more luminous clump 
stars and star No. 4 (see Fig.[Tll- 

If the observed power excesses are due to stellar oscillations 
this result shows great prospects for asteroseismology on clusters. 
However, with the unfortunate cancellation of the Eddington mis- 
sion by ESA, we might have to wait many years for a dedicated 
space project specifically aiming at asteroseismology on stellar 
clusters. Alternatively, one could initiate a ground-based network 
with high-resolution and high-sensitivity spectrographs that is able 
to detect oscillations in faint cluster stars, or a photometric multi- 
site campaign with larger telescopes (preferably larger than 2 me- 
ter, see Paper I) with stable instrumentation located at good sites 
aimed at red giants in clusters. As clearly demonstrated by our re- 
sults, the final level of the non-white noise (including drift) in the 
data will be absolutely crucial in any approach to detect solar-like 
oscillations in red giant stars. This largely favours space missions 
as atmospheric effects in ground-based data varies on time scales 
similar to the stellar oscillations. Although the Kepler and COROT 
space missions will not specifically target clusters it is expected 
that these missions will observe many red giants, some of which 
would be in clust ers. Recent results f rom the MOST satellite on the 
red giant e Oph jBarban et al. I l2007l) demonstrates that space mis- 
sions have the potential to do very interesting science on red giant 
field stars, and more K giants are on the target list for MOST in the 
future. 

We would be happy to make the data presented in this paper 
available on request. 
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Figure 10. Fourier spectra of simulated data (star No. 2) for five different 
random number seeds and two different mode lifetimes (left: 2 days, right: 
20 days). The solid white line is the smoothed spectrum. 
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